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Water famously expands upon freezing, foreshadowed by a negative coefficient of expansion of the 
liquid at temperatures close to its freezing temperature. These behaviors, and many others, reflect 
the energetic preference for local tetrahedral arrangements of water molecules and entropic effects 
that oppose it. Here, we provide theoretical analysis of mesoscopic implications of this competition, 
both equilibrium and non-equilibrium, including mediation by interfaces. With general scaling 
arguments bolstered by simulation results, and with reduced units that elucidate corresponding 
states, we derive a phase diagram for bulk and confined water and water-like materials. For water 
itself, the corresponding states cover the temperature range of 150 K to 300 K and the pressure 
range of 1 bar to 2kbar. In this regime, there are two reversible condensed phases - ice and liquid. 
Out of equilibrium, there is irreversible polyamorphism, i.e., more than one glass phase, reflecting 
dynamical arrest of coarsening ice. Temperature-time plots are derived to characterize time scales 
of the different phases and explain contrasting dynamical behaviors of different water-like systems. 



I. INTRODUCTION 

Supercooled liquids exist in a metastable equilibrium 
made possible by a separation of timescales between local 
liquid equilibration and global crystallization. pQ Super- 
cooled water is no different in this regard. However, the 
magnitude of the separation of timescales in supercooled 
water is of particular relevance due to speculation regard- 
ing the behavior of the thermodynamic properties of liq- 
uid water at very low temperatures. [2 In this work, we 
describe a theory for corresponding states that relates 
the low temperature, low pressure phase diagram with 
the time-temperature-transformation diagram for super- 
cooled water and water-like systems. Our derivations use 
scaling theories with assumptions tested against molecu- 
lar simulation. The relationships elucidate the connec- 
tions between behaviors found for different molecular 
simulation models of water and for different water-like 
substances. 

Figure [T] shows the portion of the phase diagram for su- 
percooled water relevant to this paper. Temperature, T, 
ranges from ambient conditions to deep into the super- 
cooled regime, and pressure, p, ranges from atmospheric 
conditions through the range of stability for ordinary 
hexagonal ice. Experimentally, this region corresponds 
to 150 K < T < 300 K and kbar < p < 2 kbar. [14 The 
locations of specific features relative to experiment vary 
from one molecular model to another. [15 This variabil- 
ity reflects a delicate competition between entropy and 
energy that is intrinsic to any reasonable model of water 
or water-like system. 

One manifestation of this competition is the exis- 
tence of the temperature of maximum density. We use 
T max to denote the value of this temperature at ambient 
(i.e., low pressure) conditions. For experimental water, 
T max ~ 277 K. The energy-entropy balance manifested 



in the density maximum is shifted to lower temperatures 
as elevated pressures favor denser packing. A measure 
of this shift is provided by the slope of the melting line 
or in terms of a reference pressure p Q = —AH/10AV, 
where AH and AV are, respectively, the enthalpy and 
volume changes upon melting. For experimental water, 
p « 3.7 kbar. We use T max and p Q to compare the prop- 
erties of different water models as well as to enable com- 
parison with experiment. [3H7J [14j [16] Thus, the phase 
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FIG. 1. The p-f phase diagram for supercooled water. These 
symbols refer to the pressure, p, and the temperature T, in 
units of the reference pressure, p Q , and the reference tempera- 
ture, T max , for each specific material. [3HZ] See text. The lines 
refer to theoretically expected trends as functions of pres- 
sure for the dynamical onset temperature, T Q , for the liquid 
limit of stability temperature, T s , and for a glass transition 
temperature, T g . Symbols indicate locations of these these 
temperatures in the reduced units as obtained from experi- 
ment and from different molecular simulation models, where 
the key employs standard acronyms for each system. [8 -13] For 
the glass transition, the line marks the stage where the liq- 
uid reorganization time exceeds 10 14 r o , where the reference 
time, r , is that reorganization time at the onset to correlated 
dynamics (r Q ~ 1 ps for liquid water) . 
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diagram in Fig. [T] employs the reduced variables 

P = p/Po and f = T/T max . (1) 

In this way, Fig. [I] relates results from various models 
and experiments for the onset temperature, T Q , and the 
homogeneous nucleation temperature, T s . The former, 
T Q , marks the crossover to correlated (i.e., hierarchical) 
dynamics. [T7] The latter, T s , marks the crossover to liq- 
uid instability. [18 -20J These temperatures are material 
properties. Figure [I] also shows a reduced glass transi- 
tion temperature, T g = T g /T max . This temperature is 
defined as that where the reversible structural relaxation 
time of liquid water equals 100 s. 

Glass phases of water, where aging occurs on time 
scales of 100 s or longer, are not generally accessible by 
straightforward supercooling because bulk liquid water 
spontaneously freezes into crystal ice at temperatures 
below T s . Freezing in this regime occurs in mili- second 
or shorter time scales. [21] An amorphous solid can be 
reached with a cooling trajectory that is initially fast 
enough to arrest crystallization, and finally cold enough 
to produce very slow aging. Alternatively, one may cool 
while perturbing water with surfaces that inhibit crys- 
tallization. An actual glass transition temperature of 
water is therefore not a material property because its 
value depends upon the protocol by which the material is 
driven out of equilibrium. Surface mediated approaches 
to amorphous solids can yield T g 's that are higher than 
those produced by rapid temperature quenches. The T g 
graphed in Fig. [I] is necessarily an upper bound to those 
glass transition temperatures. 

Dynamics in the vicinity of T « T s exhibits a two- 
step coarsening of the crystal phase. [20] [22] [23] First, 
disperse nano-scale domains of local crystal order form 
throughout the melt; second, on a much longer time 
scale, the nano-scale domains meld into much larger or- 
dered domains. These steps are arrested when forming 
glass. [24 Configurations appearing at the initial stages 
of this coarsening are often observed in computer simu- 
lations of water. These configurations are transient states 
that are almost as often confused with the presence of two 
distinct supercooled liquid phases, [? ] and claims that 
some water-like models do not exhibit this behavior [25] 
are based upon studies that have not examined this part 
of the phase diagram. Two distinct reversible liquids in 
coexistence would imply the existence of a low temper- 
ature critical point of the sort suggested by Stanley and 
his co-workers. [26] Not surprisingly, all reports of a low- 
temperature critical point in water or water-like systems 
locate a point on or near T s (p). In fact, all the simulation 
points clustered around that line in Fig. [I] have been in- 
correctly identified as low-temperature critical points. [8- 
[T3] Similarly, in experimental work, Mishima locates a 
putative critical point [27 close to the experimental limit 
of liquid stability. [28 We have analyzed and disproved 
this notion of a liquid-liquid transition for several differ- 
ent computer simulation models of water and water-like 
systems. [20] 



We mention the disproved notion only to emphasize 
that the equilibrium phase diagram by itself gives 
an incomplete picture of the behavior of supercooled 
water. Supercooled water, being a metastable state, 
behaves reversibly for only finite observation times. 
The specific length of that time depends on the sep- 
aration of timescales between local equilibration of 
liquid configurations and global crystallization. When 
the gap between these timescales becomes small, as it 
does for either T < T s or T « T s , time-independent 
thermodynamic properties are no longer well defined. 



II. UNIVERSAL TEMPERATURE 
DEPENDENCE OF LIQUID RELAXATION 
TIMES FOR SUPERCOOLED WATER 

In order to construct a scaling theory for the liquid 
relaxation time, rn q , we follow our previous work [19] in 
adopting a perspective of dynamic facilitation theory. [17] 
An important aspect of this perspective is that it supplies 
a universal form for the relaxation time as a function of 
temperature. This form, known as the "parabolic law", 
is 

log 10 (r liq /r o ) = J 2 (1/T - 1/T ) 2 for T < T G , (2) 

where J is an energy scale of hierarchical dynamics, T Q is 
the temperature below which that dynamics sets in, and 
t , of course, is the liquid relaxation time at the onset 
temperature T Q .[? ] This form has been used to collapse 
large and seemingly disparate collections of experimental 
and simulation data. [29] Figure [2] illustrates the nature of 
this collapse for the structural relaxation times of several 
models of water. [30H32] 

For the data shown in Fig. [2] the relaxation times have 
been calculated from the self-correlation function, 

F 8 (M) = ^e^[ ri ^- ri ^) , (3) 

for wave vectors of magnitude k ~ 27r/2.8A. Here, ri(t) 
denotes the position of a tagged molecule at time t, and 
the angle brackets stand for the equilibrium average over 
initial conditions. The time at which this F s (k,t) de- 
cays to 1/e of its initial value is defined as rii q . In all 
of the models studied, the temperature dependence of 
this time crosses over from a weak Arrhenius tempera- 
ture dependence to a super- Arrhenius temperature de- 
pendence. The location where the crossover occurs, the 
onset temperature T Q , varies from model to model as does 
the reference timescale r Q . This variability in part reflects 
quantitative differences between phase diagrams for each 
of the models. Indeed, Fig. [2}d shows that the tempera- 
ture dependence of the relaxation times can be collapsed 
by referencing the data to the temperature of maximum 
density. It is a remarkable result given the wide varia- 
tion of T max , ranging from 250 K to 320 K for the different 
models. [5] 
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FIG. 2. Low temperature variation of structural relaxation time, Tn q , for different models of water, (a) Logarithm of that time 
as a function of absolute temperature T. (b) The same data, now in units of the reference time r , and as a function of reduced 
temperature T = T/T max . The reference time, r , is the structural relaxation time at the onset temperature, T Q . (c) The same 
data, now collapsed to the parabolic form, Eq.|2] where the dashed black line is the prediction of that equation, and the dashed 
red line is the Arrhenius form that holds for temperatures above T Q . 



In addition, Fig. [2^ shows that the relaxation time data 
also collapses when referenced to the onset temperature, 
T , and that the collapsed data obeys the parabolic law 
for all T < T Q . This finding establishes that 

T Q ~ ?max • (4) 

Further, collapsing data to the parabolic law reveals that 
J/T Q varies between models of water by no more that 
5% and on average by only 1%. This universal value is 
J/T Q « 7.4. For the models considered here, r Q varies be- 
tween 0.3 ps and 8.0 ps, which largely reflects the density 
differences between models at low pressure. 

The collapse of the relaxation times for the difference 
models as a function of T max implies a universality in the 
behavior of the glass transition for these models and, by 
proxy, for experiment. As we have done previously, [19] 
we can define a locus of laboratory glass transitions as 
the locations in the phase diagram where the liquid re- 
laxation time is equal to 10 14 r o , which for many of the 
models implies TH q (T g ) « 100 s, so that with Eq. [2] we 
have, 

T g /T « (VllTo/J + l)" 1 . (5) 

Taking J/T Q « 7.4 and T Q « T max , we therefore conclude 
that for water and water-like models, T g & 0.62 T max . 
This value yields the glass-transition line graphed in 
Fig. [I] where the slope of the line is the same as that 
for T Q (p). Experimentally, the density maximum for wa- 
ter occurs at T max & 277, therefore our predicted glass 
transition is T g w 172 K. This temperature agrees with 
our previous work inferring the glass transition temper- 
ature from relaxation data of confined water. [19] It also 
provides an upper bound to values for T g obtained with 
other experimental protocols. [33] [34] 



III. MOLECULAR THEORY FOR J, r G AND T Q 

The energy, time and temperature parameters in Eq. [2] 
can be computed from microscopic theory following the 
procedures of Keys et. al.[35 The procedures are based 
upon mapping the dynamics of atomic degrees of freedom 
to dynamics of a kinetically constrained East model. [36] 
The parabolic law, Eq. [2j is a consequence of that map- 
ping. 

To illustrate the procedure for water, we have carried 
out molecular dynamics simulations of equilibrated water 
models to determine the net number of enduring displace- 
ments of length a appearing in N- molecule trajectories 
that run for observation times t ^ s . This number of dis- 
placements is 

N t ohs /At 

C a = Y.Y. e(|r,(jAt + At)-f,(jAt)|-a) (6) 

i=l j=0 

where Q(x) is 1 for x > and zero otherwise, At is the 
mean instanton time for enduring displacements of length 
a, and r^(t) is the position of particle i averaged over the 
time interval t — 5t/2 to t + St/2. The averaging over 
St coarse-grains out irrelevant vibrational motions. The 
instanton time, At, is taken to be large enough that non- 
enduring transitions are also removed from consideration. 
The two times, At > St, are determined as prescribed by 
Keys et. al.[35] 

The mean mobility (or excitation concentration) is the 
net number of enduring transitions per molecule per unit 
time, i.e., 

c a = (C a )/(Nt ohs /At). (7) 

Its dependence upon temperature and displacement 
length is illustrated in Fig. [3] According to facilitation 
theory, c a should have a Boltzmann temperature depen- 
dence, with an energy scale that grows logarithmically 
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FIG. 3. Excitation concentration for the mW model at am- 
bient pressure for enduring displacement lengthscales, a, be- 
tween 1.5 and 3.5 A. The dashed line has unit slope illus- 
trating the Boltzmann scaling, Eq. [8] The inset shows the 
logarithmic scaling of J a with a. The dashed line is a fit to 
the data for Eq.[9]with g = 0.625. 



with displacement length. That is, 

c a cxexp[-J a (l/T-l/T )] for T < T Q , (8) 

and 



J a > = J a - gJ a \n(a'/a) , 



(9) 



where a is a reference molecular length and g is a system- 
dependent constant. [? ] The data graphed in Fig. 3 
shows that for the model considered, the mW model of 
water, the theoretical expectations are obeyed. We have 
adopted the reference length a = 2.5A, which is close to 
the diameter of the molecule in the mW model, and find 
g = 0.625, T = 244 K « T max = 250 K, and J a /T Q = 23. 
According to facilitation theory [35], Eqs. [8] and [9] imply 

m(r liq /T ) = (Jl g/df) (1/T - 1/T Q ) 2 for T < T Q , 

(10) 

where df is the fractal dimension of dynamic heterogene- 



ity, which for d = 3 is about 2.6. Equation [10| therefore 
yields 



J = Ja V#/2.3df 



(11) 



where the factor of 2.3 in the square-root accounts for 
the conversion between base e and base 10 logarithms. [? 
] Applying Eq. 
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with the computed parameters yields 
J/T Q = 7.4, in good agreement with the universal empir- 
ical value reported in the previous section, that empirical 
value obtained from fitting data for various water mod- 
els. Thus, we have succeeded at deriving this value from 
a molecular calculation. 



IV. THEORY FOR CRYSTALLIZATION TIME 

To estimate the timescale for crystallization, r xt i, we 
start with the usual form, 



r^ = v-\T)e^T)/T 



(12) 



where AF(T) is the free energy cost for growing a nascent 
crystal and v~ x (T) is the timescale for adding material to 
the burgeoning phase. Typical forms for AF(T) can be 
motivated by classical nucleation theory, which has been 
shown previously to yield accurate results for nucleation 
rate of models of water at moderate supercooling. [37] In 
general, this free energy can be written as, 



AF(T)/T = $( 7 /A/i)(T/T m - 1)" 



(13) 



where 7 is surface tension for liquid-crystal coexistence, 
and $(7/ Aft) is function of the ratio of that quantity 
to the enthalpy difference between those phases. The 
ratio is approximately temperature independent. [TJ [T9] 
The temperature-dependent factor, (T/T m — l) -2 , comes 
from expanding the chemical potential difference to low- 
est non-trivial order in T — T m . 

The timescale for adding material to a growing cluster, 
v~ x (T), is expected to be relatively athermal at high tem- 
peratures, but to increase with supercooling. We expect 
^ _1 (T) ex D(T), where D is the molecular self-diffusion 
constant. Supercooled liquids generically obey a frac- 
tional Stokes-Einstein relationship, [38 



(14) 



For T > T , z = l.[39] On the other hand, for T < T G , 
z ~ 3/4. [38] This value for the exponent is predicted by 
the East model. [40] Adopting a fractional Stokes-Einstein 
relation with with Eq. [2] implies a super- Arrhenius form 
for i/" 1 ^), 



\T)(xexp 2.3 z J 2 (1/T -l/T f 



(15) 



By combining Eqs. [12}{T5] and T Q w T max , we obtain 

m(r xtl /T x ) = A (l/f - l) 2 + T (f - f m ) " 2 , (16) 

where A = 2.3z(J/T G ) 2 , T = ^>( 7 /Aft)T 2 and r G x is the 
proportionality constant in Eq. [15] 

Equation [16] can be used to fit crystallization rates in 
terms of the constants A and T. From the universal East- 
model value for z, and the universal value for J/T Q in 
water-like systems, we have A ^ 94. Further, by approx- 
imating critical nuclei as spherical and mono-disperse, 
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For the mW model we have previously determined [19] all 
the quantities on the right-hand side of Eq. [17] yielding 
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FIG. 4. Collapse of the crystallization times for models of 
supercooled water, (left) Crystallization times calculated for 
the mW model spanning the nucleation and growth regimes. 
Red markers are taken from Li et. ah [37] blue markers are 
taken from Limmer and Chandler.) 18] The dashed black line 
is given by Eq. [16] with A = 94, T = 1.3 and r* = 0.1. 
(right) Crystallization times for different models collapse. 
Black markers are our data for the mW model. Cyan and 
green markers are for two distinct sets of data for the TIP5P 
model taken from Yamada et. al.[3T] The three lines are 
Eq. [16] with A = 94 and with T and T* = xxx adjusted for 
best fits to each of the three different data sets. 



for that model T « 1.3. Equation [16] is plotted along 
side the numerical data in Fig. [4] with this parameteri- 
zation. The agreement is good over a range of 10 orders 
of magnitudes, spanning nanoseconds to seconds. The 
worst agreement is at the lowest temperature, where the 
rate is the most sensitive to the preparation of the initial 
state, as the liquid is no longer metastable at this condi- 
tion. The next section of this paper expands upon this 
point. 



V. TIME-TEMPERATURE- 
TRANSFORMATION 
DIAGRAMS 

At conditions of liquid metastability, where a free en- 
ergy barrier separates liquid and crystal basins, nucle- 
ation is the rate-determining step to form the equilib- 
rium phase. Two data sets taken from the literature have 
used different rare-event sampling techniques to com- 
pute these times for a range of temperatures for the 
mW model. Limmer and Chandler [18 computed the 
nucleation rate constant following a standard Bennett- 
Chandler procedure. [4Tj Li et. al[37] calculated the nu- 
cleation rate using forward- flux sampling [42 with an or- 
der parameter based on a crystalline cluster sizes. To 
the extent that the kinetics is nucleation limited, both 
calculations are expected to have the same temperature 
dependence. However, because each used different order 
parameters and basin definitions, the prefactors can be 
different. In order to compare both data sets, we de- 
termine the ratio of prefactors by equating the rate at 
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FIG. 5. Time scales of supercooled TIP5P water (left) and 
mW water (right) at 1 bar. Open circles are computer simula- 
tion results for structural relaxation times, rn q . Filled circles 
are computer simulation results for crystallization times, r xt i. 
One standard-deviation error estimates for the mW simula- 
tions (our results) are the size of the symbols. Error estimates 
for the TIP5P simulations |31| are unknown. The red line is 
the parabolic law, Eq. [2] and the blue and black lines are fits 
to Eq. [16] with A = 94, as predicted by theory. The grey re- 
gion is where the liquid is unstable. At the boundary of liquid 
stability, the separation between times rn q and r x ti defines the 
gap parameter, A. 



T = 220 K, which was calculated in both studies. These 
data sets are shown in left panel of Fig. [4j 

For lower temperatures, we take data sets for 
crystallization times computed from first-passage 
calculations. [43] In the calculations of Limmer and 
Chandler [18 , the first-passage time is taken from 
simulations with the mW model. Similar calculations 
by Yamada et. al[3l] based on first passage times have 
been calculated for the TIP5P model. In the latter case, 
the potential energy and structure factor were used as 
an order parameters for distinguishing crystallization. 
Both of these calculations are shown in Fig. |5j 

The crystallization times shown in Figs. [4] and [5] illus- 
trate the non-monotonic temperature dependence pre- 
dicted from Eq. 16 In the higher-temperature regime, 
nucleation rates increase because the barrier to nucle- 



ation decreases in size. In the lower-temperature regime, 
the process of crystallization is slowed by the onset of 
glassy dynamics. At conditions where the amorphous 
phase is unstable, r xt i becomes limited by mass diffu- 
sion, which from Eq. 14 is proportional to r\{ q . In this 
region of the phase diagram, the liquid state is no longer 
physically realizable. 

In plotting r x ti in Fig. [4j we use a different reduced 
temperature scale than previous plots of rn q . The dif- 
ferent scale is chosen to emphasize the crossover region, 
where nucleation and growth compete. This particular 
scale also allows for crystallization times to be collapsed 
for different models because, to first order in A/T, this 
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FIG. 6. Liquid-solid transitions of water confined to hydrophilic nanopores. (left) The phase diagram for water confined to 
hydrophilic nanopores divides into a region of liquid, glass and crystal- like states. The black line locates T m (R) while the 
blue locates T g (R). Red markers are experimental data for the melting line (circles) 48 and glass transition line (squares) . [47] 
The top part of this panel illustrates a typical configuration of ice-like water, shown in red, in contact with the hydrophilic 
nanopore, shown in grey, mediated with a premelting layer, shown in blue. The configuration is taken from molecular dynamics 
calculations of the mW model. [19] (right) The mean molecular density (top), and mean local orientational order parameter 
(bottom) , for mW water confined to a R = 20 A pore for T < T m , where g x ti is the mean value of the order parameter in the 
center of the pore. The red dashed line in the bottom right panel is the theoretical prediction from the square- gradient theory 
in Eq. [19] 



scale locates the minimum in r xt i. The location is the so- 
lution to a quadratic polynomial found by equating the 
nucleation and growth terms. [? ] This scaling holds only 
for T far below T m . Away from the singular response at 
T = T m , this form is conserved from model to model as 
it reflects the the crossover to universal structural relax- 
ation times away from the nucleation dominated regime. 

In Fig.[5]we show both rn q and r xt i to illustrate how the 
separation of timescales evolves as a function of temper- 
ature for different models of water. By considering two 
cases, where T Q w T m and where T Q < T m , we see large 
variation in the time-scale gap between liquid relaxation 
fastest crystallization. To quantify this variation between 
models, we define 



log 



10 



T x tl 



T=T S 



(18) 



and subtract Eq. [2] from [16] to predict how this gap pa- 
rameter changes with T m . For the mW model, the pa- 
rameters in this equation predict A « 1.6, in agreement 
with simulation. For experimental water, T m = 0.99, and 



T can be computed using Eq. [IT] and known values for 
7 and Aft, yielding V « 1.2. [? ] As such, Eq. 18 gives 
A = 3.4, consistent with cooling rates required to by- 
pass crystal nucleation. [21] The location of the minimum 
crystallization time for experiment can be similarly pre- 
dicted, and this yields f = 0.77 or T ^ 215 K, which is 
close to, though lower than, previous estimates. [22 One 
may also use this analysis to predict the time scales on 
which models of water will exhibit complex coarsening 
dynamics resulting in artificial polyamorphism. [20 

Table 1 summarizes materials properties noted in this 
and preceding sections. Blanks (-) in the table refer to 



properties that have not yet been determined. Viewing 
the variability between models for the values for T max and 
p elucidates how apparent different behaviors of different 
models can simply reflect different corresponding states. 



TABLE I: Summary of properties for water and wa- 
ter models, with standard acronyms identifying different 
models. [3H71 1141 1161 119] See text for meanings of symbols 
and the methods by which the properties are determined. 
Absolute temperatures are in K, pressures are in kbar and 
times are in ps. 



Model 


T 

J- max 


Po 


T m 


T 


J/T 


To T* 


r a 


Experiment 
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1.2 3.4 
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SW 
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ST2 
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TIP4P 
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0.90 
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19.4 
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VI. MESOSTRUCTURED SUPERCOOLED 
WATER 

The lengthscale over which the arguments presented in 
the previous sections are applicable to supercooled water 
reflects the lengthscale over which orientational order is 
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correlated. We have previously studied these correlations 
using the phenomenological hamiltonian of the form, 

n[q(r)} = ^(/[ 9 (r)] + ||V 9 (r)| 2 ) , (19) 

where q(r) is an order parameter that measures the 
amount of local orientational order at a point r, f(q) 
is a free energy density, 



9.1 A. This reduction in the melting temperature, 



f(q) = a(T - T s )q 2 /2 - wq 3 + uq 4 , 



(20) 



and a, w, u and m are positive constants determined by 
Aft, T m and 7. [19] This hamiltonian is isomorphic with 
that of van der Waals for liquid- vapor coexistence. [45] 
Consequently, mean profiles for q(r) subject to exter- 
nal boundary conditions yield smooth order-parameter 
profiles like those at a liquid-vapor interface. Instanta- 
neously, this field can be represented in a discrete basis 
and sampled with an interacting lattice gas. [46] Such a 
coarse-grained representation is amenable to large-scale 
computations, beyond what are tractable with atomistic 
models. 

One case of water interacting with mesoscopic inhomo- 
geneities that we have studied previously is water con- 
fined to hydrophilic nanopores. [T9] For nanopores with 
radii greater than, R > 1 nm, the properties of the water 
enclosed in the pore are sufficiently bulk-like that these 
scaling relations hold up to a perturbation due to the sur- 
face. Indeed using the expression in Eq.[5]we have shown 
that the locations of glass transitions in p Q -R plane can 
be predicted. [T9] These results are summarized in Fig. [6] 
which shows a p Q — cut through the p Q -R plane. The 
location of the the glass transition, T g for finite pores has 
been measured. [47] These points are included in Fig. [6] 
and fall on our predicted glass transition line. 

We have also previously computed the melting tem- 
perature in confinement from the partition function pre- 



scribed by Eq. 19 [19 The resulting melting temperature 
as a function of pore size and pressure is given by 

T m (p, R) = f m (p) [1 - £ m /R - £ 2 J87t(R - £ B )R] (21) 

where £ m = 27/ Aft reflects the typical spatial modu- 
lations in local order and £ s = £ m /(l — T s /T m ) is the 
renormalized length that reflects fluctuations that desta- 
bilize order. For experimental water, £ m « 2.1 A, and 



Eq. 21 , is a consequence of the silica pore wall stabilizing 
an adjacent disordered surface of water. The disordered 
surface shifts the conditions of coexistence. The melting 
line calculated from this equation is plotted in Fig. [6] and 
compared with the locations of previously determined 
freezing temperatures for water in silica nanopores. [48] 
As with the glass transition line, there is good agree- 
ment with experimental data. In out prior work, [19] we 
have also used this understanding of the phase diagram 
to explain the existence of a dynamic crossover and re- 
cent observations of hysteresis in density measurements 
for water confined to MCM-41 silica nanopores. [49] 

We mention that explanation here because it relates to 
another instance of water evolving into mesoscopic struc- 
tures, specifically the recent experimental observations 
by Murata and Tanaka. [50] Complex structure emerges 
from a mixture of water and gylcerol as it is quenched to 
low temperatures. The patterns observed depend on the 
depth of the quench and the relative concentrations of the 
two components. These patterns are reminiscent of the 
early stages of coarsening that we have found from theory 
for pure water near T s . A specific example of such struc- 
tural evolution is illustrated in Fig. [7] where the bulk free 
energy barrier to crystallization disappears. Nucleation 
occurs throughout the system and growth becomes the 
limiting timescale. This behavior is reflected in the gap 
in timescales between density and long ranged order evo- 
lution, as quantified by A. Combining the quantitative 
understanding of timescales developed in this work with 
the understanding of how ice surfaces are modulated ac- 
cording to the phenomenological hamiltonian in Eq. 19 
may admit a simple explanation for the observations of 
Murata and Tanaka. 150 
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